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Abstract 


We study the time evolution of a wave function for the spatially flat Friedmann- 
Lemaitre-Robertson-Walker universe governed by the Wheeler-DeWitt equation 
in both analytical and numerical methods. We consider a Brown-Kuchaf dust as 
a matter held in order to introduce a ’’clock” in quantum cosmology and adopt 
the Laplace-Beltrami operator-ordering. The Hamiltonian operator admits an 
inhnite number of self-adjoint extensions corresponding to a one-parameter fam¬ 
ily of boundary conditions at the origin in the minisuperspace. For any value 
of the extension parameter in the boundary condition, the evolution of a wave 
function is unitary and the classical initial singularity is avoided and replaced 
by the big bounce in the quantum system. Exact wave functions show that the 
expectation value of the spatial volume of the universe obeys the classical time 
evolution in the late time but its variance diverges. 
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1 Introduction 


Now the Inflation-Big-Bang scenario is undoubtedly a big paradigm in cosmology, sup¬ 
ported by the rapid development of observation technologies in the modern era DEI. 
Nevertheless, the singularity theorems in general relativity assert that there appears the 
initial singularity quite generically and classical physics breaks down at the very early stage 
of the universe [5]. The resolution of this initial-singularity problem requires the quantum 
description of the universe, that is, quantum cosmology. (See |1] for a review.) 

The oldest approach to quantum gravity is the canonical approach pioneered by De- 
Witt |5], based on the Arnowitt-Deser-Misner (ADM) formalism of the Einstein equa¬ 
tions |H]. Using the following ADM metric; 

ds^ = -NMt^ + hij{N^dt + dx^){N^dt + dx^), (1.1) 

where N, N^, and are functions of t{= x°) and x® {i = 1,2,3), we can write down the 
Einstein equations in the form of a constrained dynamical system. In this ADM formula¬ 
tion, the Einstein-Hilbert action Sq is written as 

Sg = j dtd^x -NH- + (total derivative). (1.2) 

Here is the momentum conjugate of hij, which is the spatial metric on a hypersurface S 
with constant t. N and act as the Lagrange multipliers corresponding to the constraints 
H = 0 and TT® = 0, respectively, where the super-momentum TT® and super-Hamiltonian H 
are functionals of hij and p®L The momentum constraints if® = 0 generate spatial diffeo- 
morphisms, while the Hamiltonian constraint H = 0 generates time reparametrizations. 

One of the possible quantizations of such a constrained dynamical system is the Dirac 
quantization [7] in which the constraint equations become operators acting on the wave 
function(al) of the spacetime d/[hij]. The resulting quantum versions of the momentum 
constraints Wd/ = 0 are formally satisfied by considering the DeWitt superspace in which 
hij take values in the quotient space of hij under the action of the group of spatial dif- 
feomorphisms. Einally, the remaining basic equation in canonical quantum gravity is the 
Hamiltonian constraint Hd> = 0, called the Wheeler-DeWitt equation. Here there is an 
ambiguity of the operator-ordering when we replace the momentum conjugates p®-^ by op¬ 
erators p®-^ := —ih6/6hij. In the present paper, we adopt the natural Laplace-Beltrami 
operator-ordering in the DeWitt superspace, originally proposed by Christodoulakis and 
Zanelli |8]. 

At a glance, the Wheeler-DeWitt equation has the form of the stationary Schrodinger 
equation with zero energy. However, since the metric in the superspace called the De¬ 
Witt supermetric has a Lorentzian (—, -|-, -f, -|-, -|-, -|-) signature at each point x® G E, the 


3 



Wheeler-DeWitt equation actually has the form of the Klein-Gordon equation in the su¬ 
perspace [S]. (See also [lU] for a review of the problem of time in quantum gravity.) This 
property of the Wheeler-DeWitt equation causes a problem of the conserved inner product. 
To avoid this problem, Brown and Kuchaf introduced ’’time” (or a ’’clock”) by matter 
helds which provide a privileged dynamical reference frame m- More precisely, they in¬ 
troduced the so-called Brown-Kuchaf dust which consists of a set of non-canonical scalar 
helds equivalent to a single timelike dust huid and cast the Wheeler-DeWitt equation into 
the form of the time-dependent Schrodinger equation. (See also [12] ■) 


In quantum cosmology, one considers only spatially homogeneous cosmological space- 
times to be quantized. In this minisuperspace approach, the metric functions depend only 
on the time coordinate, so that the resulting theory is not a quantum held theory but 
just quantum mechanics. For example, in the case of the Friedmann-Lemaitre-Robertson- 
Walker (FLRW) minisuperspace, the ehective gravitational action is 


Sg oc 


dt 



+ (total derivative). 


(1.3) 


where a = a{t) is the scale factor of the universe and p = p{t) is its conjugate. The 
initial singularity at a{t) = 0 in the classical theory is cured in quantum cosmology if the 
corresponding quantum mechanics is well-dehned. 


In [T3|, Amemiya and Koike studied the spatially hat FLRW quantum cosmology with 
a Brown-Kuchaf dust in the presence of a cosmological constant under three conceiv¬ 
able operator-orderings which are diherent from the Laplace-Beltrami one. The resulting 
Wheeler-DeWitt equation has the form of the time-dependent Schrodinger equation on the 
half-line and then the quantum cosmology is well-dehned if the Hamiltonian operator H 
acting on the wave function of the universe T = T[a] is self-adjoint or admit self-adjoint 
extension^. It was proved that, under all the operator-orderings, the Hamiltonian opera¬ 
tor admits an inhnite number of self-adjoint extensions corresponding to a one-parameter 
family of boundary conditions for the wave function at a = (@. Amemiya and Koike h- 
nally showed that the classical initial singularity is replaced by a big bounce by solving the 
Wheeler-DeWitt equation numerically with the Dirichlet or Neumann boundary condition 
in particular. 


In the present paper, we will study the same system as in |T3| but under the Laplace- 
Beltrami operator-ordering and also with a more variety of boundary conditions. In addi¬ 
tion to the study of well-dehnedness of the quantum cosmology and the initial-singularity 
avoidance, we will also clarify whether the expectation value of the spatial volume of the 
universe obeys the classical time evolution in the late time. 

^The earliest studies of the self-adjointness of the Hamiltonian operator in the self-gravitating system 
were in the context of gravitational collapse of a timelike dust shell MM- 

^The situation is similar in the system of a quantum harmonic oscillator on the half-line, in which the 
energy spectrum depends sharply on the value of the extension parameter. (See Appendix B in |16].l 
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The outline of the present paper is as follows. In section |2l the ADM formalism with a 
Brown-Kuchaf dust is reviewed. In section [31 we derive the Wheeler-DeWitt equation under 
the Laplace-Beltrami operator-ordering and determine the boundary condition for the wave 
function of the universe. Section 0] is devoted to studying the time evolution of a wave 
function. Our results are summarized in section [5l A six-parameter family of exact solutions 
to the time-dependent Schrodinger equation for a free particle or a harmonic oscillator 
obtained in im is explained in appendix A, while exact time-dependent solutions on the half 
-line constructed with the Feynman kernel are presented in appendix B. Our basic notation 
follows [18]. The convention for the Riemann curvature tensor is [Vp, 
and = R^ppy. The Minkowski metric is taken as diag(—,-1-,and Greek indices 
run over all spacetime indices. We adopt the units such that c = 1. 


2 ADM formalism with a Brown-Kuchaf dust 


We consider general relativity in the presence of a cosmological constant A in four dimen¬ 
sions, whose action is given by 

^ ^ / d^x^{R - 2A) + + Sqm. (2.1) 

where k := VSttG and G is the Newton constant. Fm is the action for matter helds and 
Sqm is the York-Gibbons-Hawking boundary term. The resulting Einstein equations are 

Rpu ~ '^dpuR + ( 2 - 2 ) 

where the energy-momentum tensor is given from S^. 


2.1 Vacuum sector 

The most general four-dimensional metric may be written in the ADM form as 

ds^ =gfj,y{x)dx^dx‘' 

= - NMt^ + hij{N^dt + dx^){N^dt + dx^), (2.3) 

where N, N^, and are functions of t{= a;°) and x^{i = 1,2,3) [6]. Let E a three- 
dimensional spacelike hypersurface with constant t and then is the induced metric on 
S. In the present paper, we assume that E is compact for simplicity. 
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In terms of the above ADM metric, the gravitational action is written as 

= J dtd^x [p^^dohij — NH — NiW) + (total derivative). (2.4) 

Here is the momentum conjugate of hij and the super-momentum if* and the super- 
Hamiltonian if are respectively given by 

if* = - 2Djp^^, (2.5) 

H (^PijP'^ - ~ 

where h := det(fij) and p := hijp^K Di and TZ are the covariant derivative and Ricci scalar 
on S, respectively. The Lapse function N and the shift vector Ni act as the Lagrange 
multipliers and, in vacuum, the Euler-Lagrange equations for N and Ni give constraint 
equations if = 0 and if* = 0, respectively. 


2.2 Matter sector 

In the present paper, we consider the Brown-Kuchar dust as a matter held m- It is a set 
of non-canonical scalar helds p, T, Z'*, and Wa {a = 1,2,3) which are equivalent to a single 


timelike dust huid, as explained below. 

The action for the Brown-Kuchaf dust is 

= d^x^pig^’^U.U, + 1), (2.7) 

where p represents the rest mass density and the one-form is dehned by 

i/^:=-(V^T) + lT,(V^Z“). (2.8) 

The Euler-Lagrange equations corresponding to T, Z'*, and Wa are 

V^pi/'*) = 0, pU^{V^Z^) = 0, V^ipWaUn = 0, (2.9) 

respectively, while the Euler-Lagrange equation corresponding to p is = —1. The 

variation 5g^'^ gives the following energy-momentum tensor; 

T^u = pU^U,. ( 2 . 10 ) 

where we have used g^'^U^Ui, = —1. From the Bianchi identity -|- Ap^**) = 0, the 

energy-momentum conservation equations = VuipU^W^) = 0 hold. 
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In terms of the ADM metric, the dynamical part of the matter action is written as 

Sm = j dtd^x (^PdoT + - AT'Pf - , (2.11) 

where P and Pa are the momentum conjugates of T and respectively m- The super¬ 
momentum Hf’ and super-Hamiltonian for the Brown-Kuchaf dust are respectively 
given by 


P°=P(V,T) + P„(V,Z“), 

(2.12) 

=^p^ + pmfHf. 

(2.13) 

The variables H4 and p are related to other ones as 


IF = - P~^P 

(2.14) 

1 =^/h\P\-^^p-‘^h^jHfHf + 1 

(2.15) 

and do not appear in the ADM form of the action (|2.11D. 



2.3 FLRW minisuperspace 

In the present paper, we will study the FLRW quantum cosmology. The line element of 
the FLRW spacetime is given by 

ds^ = —N{t)‘^dt‘^ -I- a{ty%jdx^dx\ (2-16) 

where 

(\t^ 

yjjdPdP = -- + r^(d6*^ -|- sin^ Od^p^) (2.17) 

and k = 1,0,—1 represents the spatial curvature of the universe. We have assumed the 
compactness of the spatial section of the spacetime and then the spatial volume of the 
universe is given by V{t) = Voa{t)^, where Vq := f det( 7 jj)d^a;. 

We compute 
‘S'g 

= ^ / dt!^dt{6a^N-^a) + 6Na^(J^-N-^^'^-2ANa^y (2.18) 
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where a dot denotes the derivative with respect to t. Hence, the effective vacuum action 
to give the Einstein equations for the FLRW spacetime is 


^G = 



(2.19) 


where 


Lg = 


A 1 
-w + 


A.2 


a 

—^ + 


iV2 V a2 


( 2 . 20 ) 


The momentum conjugate of the scale factor a and the super-Hamiltonian are respec¬ 
tively given by 


P 


6Lr 


6a 


QVo . 
——A aa, 




6Lg 

JW 


= ^P^ + 2VoAa^ — 6Voka 


K 


614 


Using them, we write the gravitational action in the ADM form: 


( 2 . 21 ) 

( 2 . 22 ) 


^G = 


j{pa - NH^)dt. 


(2.23) 


For the matter sector, we consistently assume = 0, H4 = 0, p = p{t), and T = T{t). 
Then we have Hf = 0 and = P, where we have taken the plus sign. Finally, the total 
action in the ADM form is given by 


S = SG + Sra= I {pa + PT- N{H^ + H^) )dt. 


(2.24) 


Now the Hamiltonian constraint = 0 is written as 

^ V + ‘AVqAo? - QVoka^ +P = 0. (2.25) 

It is noted that Amemiya and Koike considered a radiation fluid in addition in their analy¬ 
sis na. In the present paper, we don’t consider such an additional matter held for simplicity 
in order to pursue exact results as much as possible. 


The Euler-Lagrange equation for T{t) is T = 0, which is solved to give T oc t — to, 
where to is a constant. For this reason, we may use the scalar held T as a clock in the 
present system. In the comoving coordinates U^dx^ = —dt corresponding to T = t — to, 
the classical solution in this system for k = 0 with A = 0 is a(t) oc In the presence of 

positive A, the late-time behavior of the scale factor is hmi_j.oo a(t) oc We will see 

that T appears as time in the corresponding quantum system. 
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Quantum cosmology 


3.1 Laplace-Beltrami operator-ordering 

We quantize the system fl2.24|) by replacing the momentum conjugates by operators as 
p —)■ p = —ih6/6a and P ^ P = —ih{6/6T) and then the Hamiltonian constraint f|2.25p 
gives the following Wheeler-DeWitt equation: 

= (3.1) 

where 4/ = \[^[a,T] is the wave function(al) of the universe. This is the form of the 
Schrodinger equation where the scalar held T acts as a time variable. Since is an 
operator obtained from Eq. fl2.22p . there is an ambiguity of the operator-ordering. In the 
present paper, we adopt the Laplace-Beltrami operator-ordering [8] , which is natural in the 
following sense. 

By way of explanation, let us consider spatially homogeneous and anisotropic cosmolog¬ 
ical models to be quantized. According to the Bianchi classihcation, such spacetimes have 
three dynamical degrees of freedom at most, which we denote by X\t) (/ = 1,2,3). In 
this Bianchi minisuperspace, classical systems are equivalent to the dynamics of a point 
particle in a curved space. Then in general, the super-Hamiltonian fl2.22l) can be written 
in the following form: 


=^S"lX]p,pj + V[X], (3.2) 

where m is the effective mass, p/ is the momentum conjugate of , and V is the effective 
potential. From this expression, we can read off , the contravariant components of the 
supermetric Qu in the DeWitt superspace. By analogy with the canonical quantization of 
a point particle, the Laplace-Beltrami operator-ordering then requires to be 

fl° = - + 1/|A'|. (3.3) 

Here Ag is the Laplacian in the superspace: 

where Q := det(^/j). A natural inner product under this operator-ordering is 

($|T) := [ d)*T\/^d^X (3.5) 
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3.2 Wheeler-DeWitt equation 


The Laplace-Beltrami operator-ordering leads 

1 / 


= 


- 


\ 6Vo 

and finally the Wheeler-DeWitt equation fl3.1l) becomes 


1 + 21/oAa^ — GVoka^ 


6^ VqA 2 3Dnfc 


16Vo 6x'^ k'^ h? 




Jf' 


(3.6) 


(3.7) 


where x := The domain of a (and also x) is [0, oo). In terms of x, the inner product 
is simply 


($|vl/) = 


<I)*\hdx. 


(3.8) 


The equivalent Schrodinger equation to the Wheeler-DeWitt equation (13.71) is 


A ^ .AT 


where the effective potential V (x) and mass m are 

V{x 


SDo/c 2/3 

r) X X f 


m = 


8Vi 

3n2' 


(3.9) 


(3.10) 

(3.11) 


Now our problem has reduced to quantum mechanics on the half-line and we will study 
the one-dimensional Schrodinger equation fl3.9p with the inner product fl3.8l) . Hereafter we 
will consider only the spatially fiat case k = G for simplicity, in which exact solutions are 
available. 


3.3 Boundary condition for the wave function 


A quantum system is well-defined if the Hamiltonian operator 


H 


G 


2m dx'^ 


+ V{x) 


(3.12) 


is self-adjoint. Then, by the Stone’s theorem on one-parameter unitary groups, the time 
evolution of a wave function is uniquely determined by 


T(x,T) = e-*^°^/'*T(x,0). 


(3.13) 
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Let us define the Hamiltonian operator H‘^ on a dense domain T>{H^) = oo)), 

the space of smooth functions compactly supported in (0, oo), within the square-integrable 
Hilbert space £^((0, oo)). Actually, is not self-adjoint because the domain of its ad¬ 
joint operator is strictly larger than However, we may consider self-adjoint 

extensions of and then the time evolution is given by Eq. fl3.13p in their domains. 

(See [T^ [20] for a textbook and reviews.) We are going to show that our Hamiltonian 
operator with k = 0 has self-adjoint extensions characterized by one real parameter, 
namely it admits an inhnite number of self-adjoint extensions. 

For this purpose, we consider the Hilbert space eigenproblem = ±i\E' for G 

£^((0, oo)). Based on the results in Section X.l in [19], this problem reduces to solving 
the ordinary differential equations = ±i\E' for smooth square-integrable functions 

\1/. The deficiency indices n± are the dimensions of the space of general solutions to these 
differential equations without imposing boundary conditions at a: = 0. Namely, n_|_ and 
? 7 ,_ denote the numbers of parameters contained in the general solutions to the differential 
equations with +i and —i, respectively. If n± = 0, then is essentially self-adjoint and 
no further boundary conditions are required. If n+ ^ ri-, then has no self-adjoint 
extensions and the quantum system is ill-dehned. If n_^ = n_ ^ 0, then has self-adjoint 
extensions which require the imposition of further boundary conditions with n+ = n_ 
parameters. 

In the case of A = fc = 0, the general solutions to the differential equations 
are 

T/ \ ^ f I l±i \ „ f I 16Ho l±i \ m i a\ 

Independent of the sign in the differential equation, the term with Ci diverges for a; —)■ oo 
and only the term with C 2 provides a solution in £^((0, 00 )). This shows that the deficiency 
indices are n+ = ?t._ = 1 and therefore our Hamiltonian operator with A = A; = 0 has 
an inhnite number of self-adjoint extensions characterized by one parameter. 

The result is the same also with a positive cosmological constant (A > 0). In this case, by 
the scaling transformations z := /da: and T := (3K^h/d^/16Vo)T with/d := (64V))^A/3/p^h^)^/^, 
the Schrodinger equation fl3.9p becomes iLT = id'^/dT, where 

(3.15) 

Hence, we consider the following ordinary differential equations 

(3.16) 

for smooth square-integrable functions T {z) without imposing boundary conditions ai z = 0 
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(and hence x = 0). The general solutions to these differential equations are 

f (2) = H) + /_,/4(±fc74)Vc4lA3/4(±'2V4) - Ky,{±iz‘l4,)\ j, 

(3.17) 

where Cz and are constants and /^(x) and Ka{x) are the modified Bessel functions of 
the first and second kinds, respectively. (See Section 12.14 in [21].) Clearly, the term with 
the constant C 3 diverges for z ^ 00 and only the term with the constant C 4 provides a 
solution in /1^((0, cxo)). 


We have shown that our Hamiltonian operator with k = 0 has self-adjoint exten¬ 
sions which require a boundary condition with one parameter. This boundary condition 
is determined so as to satisfy the symmetric property = H^. Using the Schrodinger 
equation fl3.9p and integration by parts together with the fall-off condition at infinity, we 
obtain 




2 m \ dx 



a :=0 


(3.18) 


In order for the surface term to be vanishing, we impose boundary conditions on 4) and T 
such as 


$(0,T) + L—(0,T) = 0, 
T(0,T) + l|^(0,T) = 0, 


(3.19) 

(3.20) 


where L is a real constant, and then is realized. These boundary conditions 

also emerge from the deficiency space machinery as discussed in [T9] . 


Since the boundary condition fl3.20l) contains one real parameter L, the Hamiltonian 
operator admits an infinite number of self-adjoint extensions and each value of L gives 
a different quantum system. This boundary condition ensures unitarity c}(T|T)/9T = 0, 
shown as 


A 


(>i>i>i>> 


2 m \ dx dx J 


= 0 . 


x=0 


(3.21) 


L = 0 and L = 00 correspond to the Dirichlet and Neumann boundary conditions at x = 0, 
respectively, and other values of L correspond to the Robin boundary condition. 


4 Unitary evolution of the quantum universe 


In this section, we will solve the the Schrodinger equation (13. 9 p (which is equivalent to 
the Wheeler-DeWitt equation (13.7p ) with k = 0 and see dynamical properties of a wave 
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function depending on the extension parameter L and A. In particular, we will check 
whether the expectation value of the spatial volume of the universe, which is proportional to 
(a^)(= (x^)), obeys the classical time evolution (a?) oc (for A = 0) or (a^) oc exp(\/M^T) 
(for A > 0) in the late time. 


In the case of the full-line x G (—cxo, cxd), the following Ehrenfest’s theorem holds; 

/ dV 

dT"^ \ dx 


m- 


(4.1) 


and therefore (x) follows classical orbits. The above equation is modihed in the half-line 
case X G [0, oo). Using the Schrodinger equation (13.911 and integration by parts, we obtain 

d‘^{x‘^) h d 


m- 


dT^ 


2idT 

qh^ 
+ — 
4m 


x^ 


X 


-T — 4/ — 

dx / j 0 

^ ‘ -'h - 2 — -^ 1 -(g - 


dx'^ 

°° dV 

qx<i 


dx"^ 


dx dx 


d^\ 

dx dx J 


dx 


q{q — l)h^ 


POO 


4m 


{q-2)x‘^ —- + 4x'> Idx 


dx 


dx dx 


(4.2) 


for constant q. The surface term at a; = 0 in the hrst term vanishes for positive q by the 
boundary condition fl3.20p . Then, under the assumption that the surface terms at inhnity 
vanish, Eq. (I4.2p reduces to 


m- 


d'^{x'^) 

dT^ 


qh^ 

4m 


x^ 


d‘^4> 

“ ' Vl> + Vl/*-2- 


dx"^ 

°° dV 

qx^ i\]>*__\]>da; 


dx"^ 


dx dx 


— {q — l)x^ ^ 


T -F T*— I 


dx 


dx J 


x=0 


dx 


q{q — l)h^ 
4m 


POO 


{q-2)x‘’ - - + 4x‘^ Idx 


for positive q. For q = 1, this becomes 


m- 


(x) 

dT^ 




- 


Am \ dx'^ 


T T 


dx 

dx'^ 


- 2 - 


dx dx 

aT* dAl 


dx dx 


x=0 


dV 

dx 


(4.3) 


(4.4) 


and therefore in general, the expectation value of x does not follow classical orbits. 


Our interest is the evolution of (a^)(= {x"^)) and this quantity satishes the following 
equation: 


d^x^) r rd^Af* ^^d^AJ dAJ*dAi 

5 T2 2m [ \dx^ dx^ dx dx 

2!fdA/*dA/ 

-o- 

m dx dx dx J 


dA!* 

-4> -\- v]/*- 

dx dx 


x=0 


(4.5) 
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4.1 Analytical results 


4.1.1 Exact wave function I 


As seen in Eq. fl3.9p . the Wheeler-DeWitt eqnation fl3.7p with A > 0 and fc = 0 is eqniv- 
alent to the Schrodinger equation for a harmonic oscillator with negative mass. By the 
transformation x = (8Vo)x, Eq. fl3.9p with k = 0 becomes 


I72-2T 


(4.6) 


where k"^ := 3A/4. In the case of the Dirichlet or Neumann boundary condition at the 
origin, exact time-dependent solutions to Eq. fl4.6p are available. 

First, we use the six-parameter family of exact solutions obtained in |T^. (See Ap¬ 
pendix ISl) The solution to Eq. (I4.6p is obtained by x —)■ x and t —T from Eqs. flA.2l) 
and flA.l 6 p - flA. 22 p as 


T(x,T) =T„(x,T) 

gj(a(T)S^+5(r)x-|-K(T))+i(2n+l)7(T) 


g-(/3(T)x+£(T))V2jy^ (/?(T)x + eiT )), (4.7) 


^2"n!/i(T)v^ 

where Hn{x) (n = 0,1, 2, • • •) is the Hermite polynomials and 

/i(r) =^iQ\J /Sq sinh^ kT/kP' -|- (2q;o sinh kT — cosh fcT)^, 
kao cosh.2kT — (sinh2A;T/fc)(/dQ — Ak'^al + fc^)/4 
/Sq sinh^ kT/k? -[- (2q;o sinh kT — cosh kT^ 


a{T) 

m 


/9q sinh^ kT/k'^ + {2ao sinh kT — cosh kTy 


1 f P^sinh kT/k 

7(1 I = 7 o- arctan - =-= — 

2 V 2tto sinh kT — cosh kT 


6 {T) = - 
e{T) = - 

k{T) =Kq + 


5o(2Q;osinh kT — cosh kT) -f- sinh kT/k 
/3q sinh^ kT/k'^ + (2q;o sinh kT — cosh kT)‘^ ’ 

£o(2do sinh kT — cosh kT) — /3 q6o sinh kT/k 

sinh^ kT/k'^ + {2ao sinh kT — cosh kT)^ 

sinh^ kT e^f^fkaoeo — fo6o) — 

/do sinh^ kT/k'^ + (2q;o sinh kT — cosh kT)'^ 
1 sinh 2kT — 6 q 

4 k /do sinh^ -I-(2Q;osinhfcT — coshfcT)^ 


(4.8) 

(4.9) 

(4.10) 

(4.11) 

(4.12) 

(4.13) 


(4.14) 
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One of the seven parameters po) ctO) f^Oi 7o, niay be used for normalization, so that 

the number of independent parameters is six. In the free-particle limit (fc —?• 0), we have 


/i(T) + 1, 

/^O 


m 

6 {T) 


V/So^-^ + l’ 

<^o ~ ^oI^qT 
'' /3o^T2 + 1 ’ 



where the parameter ao has disappeared. 


^ ’ 2(/d4T2 + l)’ 

(4.15) 

7 iT) = 7o + ^ arctan(^oT), 

(4.16) 

^0 + (do^oT 

^ y//3o^T2 + l’ 

(4.17) 

- QT 

2 (/do^T2 + 1) ’ 

(4.18) 


This family of solutions with £o = <5o = 0 (hence e{T) = 0) are also solutions in cxd)) 

with the Dirichlet or Neumann boundary condition at the origin. The Dirichlet boundary 
condition is satished for odd n, while the Neumann boundary condition is satished for even 
n. The solutions are regular everywhere for T G (— cxo, cxd) and hence the classical initial 
singularity is avoided in the corresponding quantum system. Figure [1] shows the evolution 
of ITF for n = 1 and n = 2. 


Let us discuss the late-time evolution of the expectation value of the spatial volume of 
the universe, which is proportional to (a^)(= For e(T) = 0 {eq = 6 o = 0) and an 

integer q, we obtain 

1 

/ T)|^da; 


2”n!/i(T)v^7o 

1 


^qg-/3(T) X Hn{(3{T)xfdx 

pzLoo 


y^e y Hn{y) dy, (4.19) 


2-n!/i(T)/3(T)i+^0F7o 

where the sign in ±cxo corresponds to the sign of /do. Using this, we compute 
xy\d/n{x,T)\‘^dx _ 1 


r‘diOD 


(x'^) oc 


where we used 


\dI^{x,T)\^dx 2-^n\/3{Ty^ J, 


y^e-y H^{yYdy, (4.20) 


e y^HYyfdy = v^2’" ^n\. 


(4.21) 


Equation (I4.2np shows that the time-dependence of (a^) is given by 

r +1) 

/3(T)2 I /d(C^ ( sinh^ kT -|- (2q;o sinh kT — cosh kT^ 


(a^) oc 


for A = 0, 
for A > 0. 


(4.22) 
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(al) (a2) 




Figure 1: Time evolution of |\I/p for the exact wave function fl4.7p with (al) A = 0 and 
n = 1, (a2) A = 0 and n = 2, (bl) A = 0.5 and n = 1, (b2) A = 0.5 and n = 2, where we 
set /io = 1, do = Po = 1, 7o = 1; = 0, Kq = 1, and Eq = 0. The big bounce occurs at 
T = T-q and |\I/p is time-symmetric with respect to Tb, where Tb is dehned by Eq. 04.231) . 
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In both cases, (a^) is positive definite and it admits only one local minimum at T = Tb, 
where 


Tb = 


0 

1 + ^^(1 + 2do)^A 

4l V/9o+ fc2(l-2do)V 


for A = 0, 
for A > 0, 


(4.23) 


corresponding to the transition time from the contracting phase (d(a^)/dT < 0) to the 
expanding phase (d(a^)/dT > 0). Thus, the classical initial singularity is avoided and 
replaced by the big bounce at T = Tb in the quantum system. Figure [2] shows the function 
/?(T)-2((x (a^)). 


1/ p(Th 



Figure 2: The function /9(T)“^(oc (a^)) for A = 0 (thin curve) and A = 0.5 (thick curve), 
where we set do = 1 and /do = 1. The big bounce occurs at T = Tb, defined by Eq. (I4.23p . 
where Tb = 0 for A = 0 and Tb — 0.47253 for A = 0.5. 


In addition, the late-time evolution of (a^) is 


. y . T^ for A = 0, 
hm (a ) o(. < nur r ^ 

T-^oo' I for A > 0. 


(4.24) 


Since we have k = ^3A/4 in our system, this proves the convergence to the classical 
evolution in both cases. On the other hand, the variance of diverges as 


lim Var(a^) = lim ( (a^’) — (a^)^ ) oc hm /d(T) ^ cxd. 

T—>-oo T^ooV / T—^oo 


(4.25) 
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4.1.2 Exact wave function II 


To the Schrodinger equation (I3.9p with A > 0 and /c = 0, there is also another class of 
exact solutions with Dirichlet or Neumann boundary condition, constructed by using the 
Feynman kernel. (See Appendix [B] for derivation.) 

In the case of £^((—cxo, oo)), the typical initial profile is the following Gaussian wave 
packet: 


T(3:,0) = Gexp^^- ^^ ^ 

where Pq, Xq, and a are constants and the normalization constant C is l/(27rcr^)^/^. This 
wave function represents a moving wave packet with its peak at a: = a:o and its momentum 
p, as shown by (a;)(0) = Xq and (p)(0) = po- 


In the half-line case £^([0, oo)), the physical meanings of the constants Xq and po are less 
clear and the profile (I4.26p does not satisfy neither the Dirichlet nor Neumann boundary 
condition at a: = 0. From this observation, we modify the initial profile as 


T(a;, 0) = xexp 


f {x- Xq)^ .Pqx] 

\ ^ h ) 


for L = 0 and 


T(a;,0) 




[ {x-xof .pox\ 


(4.27) 


(4.28) 


for A 7 ^ 0, which satisfy the boundary condition fl3.20p . 


The solution with the Dirichlet boundary condition and the initial prohle fl4.27p is ob¬ 
tained from Eq. flB.12p by the reparametrization u ^ iu) and identifying m = —8Vo/3k^ 
and = 3A/4: 


T(x,T) = — 


mu 


2p V ^nihsinhuT 
1 
2 


exp 


Xn imux"^ cosh uT \ 


4(j^ 


+ 


2fisinh uT 


■ 


where F{z) is the Dawson function defined by 




2^^ 


(4.29) 


F{z) :=e-"' [\^"dw 

Jo 


(4.30) 
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and rj = rj{T), ( = ({x,T), and ( = C{x,T) are given by 

1 imoj cosh uT 

r]{T) 


4(j2 2hsmh. ujT 
Xo , .Po imux 

C(X,T) - J-r^, 

y/ Xo .Po imux 

+ jr-^. 


(4.31) 


On the other hand, the solntion with the Nenmann bonndary condition and the initial 
prohle fl4.28p (with L —)■ cxd) is obtained from Eq. flB.lSp as 


4/(x,r) = 


mu 


■ exp 


X 


xl imux"^ cosh uT \ 

27rfhsinhci)T ^ 4 ct^ 2hsinha)T / 

J_ / 1 , [^f (^gCV4i? ^gCV4^?^ . -^T^f 


2 p\2y 7] 

Xo , -Po 


+zCF 


2v^ 


+iCF 


K 


-1 


,CV4i7 


_j_ g 


CV4i7^ 




iC, \ ^ 


+2 


+ - 


2^/h/ VV 


K 


(4.32) 


In the limit to the zero cosmological constant A —)■ 0 (a; ^ 0), the solntions (I4.29p and 
fl4.32p rednce to 


^{x,T) = 


m 


^—XQ/4:a^-\-imx^/2hT 


2p Y 2'KihT 
1 Fk 


X 




Cf(- 


V 2^1) 


and 

^(a:,T) = 


^—XqI -\-i7nx^ l2hT 


2tt ihT 


X 




■ 4 ■ 


2 p\ 2 y 7] 

Xo , .Po 
+ 7 


+tCF 


K 


2v^ 




2(j^ h 


-1 


1 / 77 ( /-2/yi^ ^2 j \ 


2 \j p 


gCV4. + gC- 


+ - 


(-^)+2 

V 

\+—F 


respectively, where complex fnnctions p, (, and ( are now 

1 im 


V(.T) 


4(t2 2hT' 


M T\ ^0 , Po '^onx 
C(x.T) 

C(x.T) 


(4.33) 




J y/v \ ‘^y/nj VF V 2^ 

(4.34) 


(4.35) 
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Although it is difficult to see analytically the bouncing behavior of (x^)(T)(= (a^)(T)) 
in these solutions, we can evaluate its asymptotic behavior for T oo. For the solution 
with the Dirichlet boundary condition fl4.29p . changing the coordinate as s := ux/ sinhcuT, 
we obtain 


T)|^da; 


m fsinhuTV 


8 nh\ri\‘'^ 


CO 






'^Vv 


-ICF 

-iCF 


K 


2^/h 

K 

' 2 ^ 


where 

1 imo) cosh wT xq po im - xq po im 

”= 4 ^- 2fismhi,r ' ^ = ^ + = ^ + 

Since 14^(0:, T)pda: is constant, we evaluate 


ds, 

(4.36) 

(4.37) 


lim 

T^oo 


x^|\['(x, T)pdx /sinhco’TX'^ 

jyi>i.(..r)pd. 


(4.38) 


under the assumption that the limit commutes with the integral. We obtain the same 
result for the solution with the Neumann boundary condition (I4.32p . In the limit of A —)■ 0 
(a) ^ 0), we obtain 


in both cases. 


lim (x'^} oc T'^ 

T^oo 


(4.39) 


Therefore, the late-time evolution of (a^) is 


lim (a^) oc 

T^oo 


for A = 0, 
for A>0. 


(4.40) 


Since we have to = 1^3A/4 in our system, this again shows the convergence to the clas¬ 
sical evolution. Unfortunately, the Feynman kernel is not available to construct solutions 
with the Robin boundary condition. We will study such solutions numerically in the next 
subsection. 


4.2 Numerical results 

In the previous subsection, we studied exact solutions satisfying the Dirichlet or Neumann 
boundary condition at the origin. However, dynamical properties of solutions with the 
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Robin boundary condition are still not clear. Here we study this problem by solving the 
Schrodinger equation fl3.9p with k = Q numerically from the initial prohle fl4.27p or fl4.28p . 
We treat L as a parameter controlling the boundary condition with hxed values of a, Xq, 
and po and adopt the Planck unit G = h = c = 1. In addition, we set Vq = 1 which means 
that the spatial volume of the universe is the Planck volume when a (and hence x) is unity. 

In our numerical calculations, we set the space step Ax = 0.05 and time step AT = 
0.05^/4. We conhrmed that even with the smaller stepsize the results in the hgures and 
tables are unchanged. We also conhrmed that the exact solution fl4.7p with n = 1 and 2 
can be constructed numerically in the period of time shown in Fig. [1] In the long-time 
computations, however, we could not keep enough accuracy to verify the convergence of 
{a?) to the classical evolution. 

We solved the Schrodinger equation (I3.9p numerically by Maple with xq = 20, po = 0.8, 
and (7 = 5. The time evolutions of |\kp(x,T) with A = 0, 0.1, 0.5 and 1.0 are shown in 
FigsISHSl respectively. The wave function is regular everywhere during the evolution and 
it contracts initially and moves back after some moment. The big bounce behaviors of 
(a^)(= (x^)) are shown in Tables [THll For a hxed value of A, the diherence in the prohles 
of iTp with diherent values of L appears only around the origin. As A increases, the big 
bounce occurs sooner and there appears less oscillation around then. 


Table 1: Values of (x^)(= {a?)) for various values of L in the case of A = 0 with xq = 20, 
cr = 5, and po = 0.8. 


T 

T = 0 

L = ±oo 

T = 10 

L = -10 

T = 20 

L = -20@ 

0 

522.06 

524.08 

524.41 

523.62 

524.26 

523.86 

0.5 

368.02 

370.06 

370.30 

369.72 

370.19 

369.90 

1.0 

242.89 

244.78 

244.94 

244.53 

244.87 

244.67 

1.5 

146.66 

148.23 

148.33 

148.06 

148.29 

148.15 

2.0 

79.332 

80.407 

80.550 

80.232 

80.482 

80.323 

2.5 

40.910 

41.319 

41.777 

40.880 

41.547 

41.096 

3.0 

31.392 

30.963 

32.273 

29.748 

31.612 

30.337 

3.5 

50.779 

49.338 

52.123 

46.729 

50.724 

47.997 

4.0 

99.069 

96.446 

101.19 

91.949 

98.814 

94.140 

4.5 

176.26 

172.29 

179.26 

165.61 

175.78 

168.87 

5.0 

282.36 

276.86 

286.19 

267.86 

281.54 

272.26 

5.5 

417.36 

410.16 

421.92 

398.78 

416.07 

404.35 
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(a) |XJ/|2 |XJ/|2 



Figure 3: Time evolution of |Tp before (left) and after (right) the big bounce (around 
T = 3.0) for (a) L = 0, (b) L = ±oo, (c) L = 10, and (d) L = —10 with A = 0 and Xq = 20, 
cr = 5, and po = 0-8. 
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Figure 4: Time evolution of |\l/p before (left) and after (right) the big bounce (around 
T = 3.2) for (a) L = 0, (b) L = ±oo, (c) L = 10, and (d) L = —10 with A = 0.1 and 
Xo = 20, a = 5, and po = 0.8. 
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Figure 5: Time evolution of |Tp around the big bounce (around T = 1.2) for (a) L = 0, 
(b) L = ±cxo, (c) L = 10, and (d) L = —10 with A = 0.5 and Xq = 20, ct = 5, and po = 0.8. 
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Figure 6: Time evolution of |\I/paround the big bounce (around T = 0.6) for (a) L = 0, (b) 
L = ±cxo, (c) L = 10, and (d) L = —10 with A = 1.0 and xq = 20, a = 5, and po = 0.8. 
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Table 2: Values of (x^)(= (a^)) for various values of L in the case of A = 0.1 with xq = 20, 
cr = 5, and po = 0.8. 


T 

L = 0 

L = ±oo 

L = 10 

L = -10 

L = 20 

L = -20@ 

0 

522.06 

524.08 

524.41 

523.62 

524.26 

523.86 

0.5 

375.84 

377.92 

378.16 

378.06 

377.57 

377.76 

1.0 

267.37 

269.42 

269.59 

269.52 

269.15 

269.30 

1.5 

188.46 

190.38 

190.52 

190.46 

190.17 

190.29 

2.0 

133.17 

134.86 

135.01 

134.94 

134.65 

134.76 

2.5 

97.318 

98.658 

98.991 

98.828 

98.304 

98.483 

3.0 

78.206 

79.045 

79.877 

79.458 

78.263 

78.645 

3.5 

74.390 

74.543 

76.325 

75.424 

72.923 

73.703 

4.0 

85.582 

84.810 

88.095 

86.435 

81.853 

83.270 

4.5 

112.63 

110.62 

116.05 

113.31 

105.75 

108.08 

5.0 

157.56 

153.92 

162.24 

158.05 

146.45 

150.02 

5.5 

223.77 

217.97 

230.11 

224.00 

207.09 

212.28 


5 Summary 


In the present paper, we have studied the time evolution of a wave function for the spatially 
flat FLRW universe governed by the Wheeler-DeWitt equation, with or without a positive 
cosmological constant A. We have adopted the Laplace-Beltrami operator-ordering and 
considered a Brown-Kuchaf dust as a matter held. Then the system has reduced to quantum 
mechanics on the half-line and the Wheeler-DeWitt equation has the form of the time- 
dependent Schrodinger equation for a harmonic oscillator with negative mass, where a 
scalar held T acts as a time variable. 

Self-adjoint extension of the Hamiltonian operator admits a one-parameter family of 
boundary conditions at the origin in the minisuperspace. For any value of the extension 
parameter L, the time evolution of a wave function is unitary and the corresponding quan¬ 
tum system is totally well-dehned. We have shown that the classical initial singularity is 
avoided and replaced by the big bounce in the quantum system. These properties have 
been shown also under the diherent operator-orderings in |13j . 

We have also studied the problem whether the expectation value of the spatial volume of 
the universe (a^) obeys the classical evolution in the late time. This is a nontrivial problem 
because the Ehrenfest’s theorem is not valid in quantum mechanics on the half line. We 
have used exact solutions with the Dirichlet or Neumann boundary condition at the origin 
and analytically showed the convergence to the classical evolution both in the cases with 
and without A. However, this does not mean the classicalization of the quantum universe 
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Table 3: Values of (x^)(= (a^)) for various values of L in the case of A = 0.5 with xq = 20, 
cr = 5, and po = 0.8. 


T 

L = 0 

L = ±oo 

L = 10 

L = -10 

L = 20 

L = -20@ 

0 

522.06 

524.08 

524.41 

523.62 

524.26 

523.86 

0.2 

464.18 

466.26 

466.55 

466.42 

465.84 

466.06 

0.4 

423.19 

425.36 

425.63 

425.51 

424.97 

425.18 

0.6 

396.62 

398.93 

399.19 

399.08 

398.56 

398.76 

0.8 

382.88 

385.38 

385.63 

385.52 

385.01 

385.21 

1.0 

381.13 

383.89 

384.14 

384.03 

383.51 

383.72 

1.2 

391.27 

394.37 

394.62 

394.51 

393.97 

394.19 

1.4 

413.92 

417.44 

417.72 

417.60 

417.01 

417.25 

1.6 

450.43 

454.50 

454.81 

454.67 

454.02 

454.28 

1.8 

502.99 

507.78 

508.12 

507.97 

507.23 

507.53 

2.0 

574.77 

580.47 

580.87 

580.70 

579.83 

580.18 

2.2 

670.08 

676.94 

677.42 

677.22 

676.18 

676.60 


in the late time because the variance of is diverging. 

Unfortunately, exact solutions are not available in the case of the Robin boundary con¬ 
dition. In order to verify the convergence in such cases, long-time numerical computations 
with high accuracy are required. Such numerical studies will also clarify the effect of the 
spatial curvature of the universe on the evolution. These are left for future investigations. 

Another promising direction of future research is to generalize our study in a more gen¬ 
eral minisuperspace such as the Bianchi minisuperspace. In such cases, the corresponding 
quantum systems are higher-dimensional and the problems of the self-adjoint extension and 
initial-singularity avoidance are highly nontrivial. Also, quantization of the inhomogeneous 
universe is another possible direction. Classically, the general spherically symmetric solu¬ 
tion with a dust fluid is the Lemaitre-Tolman-Bondi solution. The corresponding quantum 
system is inhnite dimensional, namely a quantum held theory, and therefore a totally dif¬ 
ferent treatment of the system is needed [22]. Those results could shed light on generic 
properties of canonical quantum cosmology and will be reported elsewhere. 
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Table 4: Values of (x^)(= (a^)) for various values of L in the case of A = 1.0 with xq = 20, 
cr = 5, and po = 0.8. 


T 

L = 0 

L = ±oo 

L = 10 

L = -10 

L = 20 

L = -20@ 

0 

522.06 

524.08 

524.41 

523.62 

524.26 

523.86 

0.1 

492.70 

494.75 

495.06 

494.31 

494.92 

494.55 

0.2 

471.45 

473.56 

473.87 

473.14 

473.73 

473.37 

0.3 

457.70 

459.89 

460.19 

459.48 

460.05 

459.70 

0.4 

451.02 

453.32 

453.61 

452.91 

453.48 

453.13 

0.5 

451.20 

453.66 

453.95 

453.24 

453.82 

453.46 

0.6 

458.27 

460.91 

461.20 

460.48 

461.07 

460.71 

0.7 

472.43 

475.29 

475.59 

474.84 

475.46 

475.08 

0.8 

494.10 

497.23 

497.55 

496.75 

497.41 

497.01 

0.9 

523.93 

527.39 

527.73 

526.88 

527.58 

527.16 

1.0 

562.82 

566.67 

567.04 

566.12 

566.87 

566.42 

1.1 

611.92 

616.25 

616.64 

615.64 

616.47 

615.97 
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A Six-parameter family of exact solutions 


For the following Schrodinger equation for a harmonic oscillator with positive mass; 

1,2 2 . 

- , 

2 dx^ 2 dt 


(A.l) 


there is a six-parameter family of exact solutions [17)4 : 






g-(/3(qx+£(t)) /^HMt)x + e{t)), (A.2) 


^We have done coordinate transformations and reparametrization from the original expressions in 













where Hn{x) is the Hermite polynomials and 


(A.3) 

(A.4) 


/i(t) 

a{t) 

m 

lit) 

m 

e{t) 

K{t) 


=/io Y /Sq kt/kl + (2q;o sin kt + cos kty, 

kao cos2kt + {sm2kt/k){f3Q + Ak'^a'^ — k‘^)/A 
PqSiv? kt/k"^ + {2aQ sin kt + cos kty ’ 

__^0_ 

sin^ kt/kl + (2q;o sin kt + cos kty 

1 / jl/sinkt/k \ 

2 \ 2q;o sin kt + cos kt J 
So{2ao sin kt + cos kt) + eo/^o sin kt/k 
13Q srn^ kt/k"^ + {/2aQ sin kt + cos kt)'^ ’ 

So (2 q;o sin kt + cos kt) — //qSq sin kt/k 
sj/3 q sin^ kt/k'^ + (2q;o sin kt + cos kt)'^ 

sin^ kt £ 0/^0 ~ /^o^o) ~ kaoE^ 

° kl //^sin^ kt/k"^ + {/2ao sin kt + cos kt)"^ 
lsin2/ct Sq/Sq — (5o 


4 k (3^ sin^ kt/kl + {2aQ sin kt + cos kt)"^ 


(A.5) 

(A.6) 

(A.7) 

(A.8) 

(A.9) 


One of the seven parameters /io, cto, /do, 7o, <^o, ^ 0 , may be used for normalization, so that 
the number of independent parameters is six. The norm of in the case of the full line 
is given by 


where we used 


\^nix,t)y<3x 


1 

ho/do 


kkjnix) Hnix)c 


(3x = ^/^^2"'n\6„ 


(A. 10) 


(A.ll) 


The solution for a free particle is obtained in the limit /c —)■ 0 as 

t/3t/2 

/d^f^ + 1 ’ 


/^(^) =/aoV/^o^^ + 1, a{t) = 


m 


m 


/do 1 - 

^^====, 7(^) = 7o - X aTctan{f3^t), 

VPot +1 ^ 


<^0 + £o/3ot 

' P^ot^ + 1 ’ 


e{t) = 


£0 — /3oSot 

V/do^t' + l’ 


,, ,2 ^o/3lSo 1 el/3^-6^ 

■ ‘ MfiTT^ 7WTT’ 


(A.12) 

(A,13) 

(A.14) 

(A.15) 
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where the parameter oq has disappeared. By the reparametrization k = ik and ia^ = Oq 
in Eqs. flA.3p - flA.9p . we obtain the solution for an inverted-harmonic oscillator as 


a{t) 

m 


=/io \J /Sq sinh^ kt/k"^ + (2ao sinh kt + cosh kty, 
kaQ cosh2fct -|- (sinh2H/fc)(/dQ — Ak'^al -|- fc^)/4 
/5 q sinh^ kt/k^ + {2ao sinh kt + cosh kty 
_ _ ^0 _ 

\J /5q sinh^ kt/k'^ + (2q;o sinh kt + cosh ktY 


7(^) 

m 

e{t) 

K{t) 


1 / I3nsmh.kt/k \ 

7o -arctan -, 

2 \ 2ao sinh kt + cosh kt J 

5o(2do sinh kt + cosh kt) + sinh kt/k 
/3q sinh^ kt/k^ + {2ao sinh kt + cosh kt) ^ ’ 

£o(2do sinh kt + cosh kt) — /Sq^o sinh kt/k 

/5q sinh^ kt/k"^ + (2q;o sinh kt + cosh kt)"^ 

sinh^ kt eo/do(fcdo£o — PoSq) — fcdo^o 
° k"^ /dg sinh^ + (2Q;osinh fct + coshfct)2 

1 sinh 2kt ^q/Sq — 

4 k sinh^ + (2Q:osinh H + cosh 


(A.16) 

(A.17) 

(A.18) 

(A.19) 

(A.20) 

(A.21) 


(A.22) 


The solution of the Schrodinger equation fl4.6l) for a harmonic oscillator with negative mass 
is obtained from the above by the transformation t —)■ —T. 


B Exact solutions constructed with the Feynman ker¬ 
nel 


In this appendix, we present exact time-dependent wave functions on the half-line con¬ 
structed with the Feynman kernel. In the case of the full line, the time-dependent solution 
of the Schrodinger equation 


2m dx^ 


+ I7(x)$ = ih 




with an initial prohle <h(x, 0) is given by 


<h(a;, t) 


K{x, t] x', 0)<I)(a;', 0)da;'. 


(B.l) 


(B.2) 
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For a free particle or harmonic oscillator, the Feynman kernel K{x,t]x',Qi) is known |23j . 
In the case of the harmonic oscillator V{x) = {l/2)mw‘^x‘^, it is given by 


K{x,t;x',0) 


mu fimulix"^ + x'"^) cos ut — 2xx'} 

exp - 


27iihsmut \ 2hsmut 

Taking the limit a; ^ 0, we obtain the Feynman kernel for a free particle: 


(B.3) 


, , I m f imix — x'Y 


which is also written as 


1 

K(x,t;x',0) =— / 
27t 


(B.4) 


(B.5) 


Using them, we can construct time-dependent solutions on the half-line with the Dirichlet 
or Neumann boundary condition at a: = 0. The solution is then given by 


^{x,t) = / K{x,t-,x',0)^{x',0)dx' 


(B,6) 


and the expressions of K{x,t]x',0) will be shown below. For derivation, we will use the 
following Formula 7.7.3 in 


+ ^F{z/^) (Reiv) > 0 ), 

2\ T] ^ 


(B.7) 


where F{z) is the Dawson function dehned by 


F{z) := e-^ / e"' dw. 


(B.8) 


We will also use the following: 


X e 


'0 


°° 9 1 rl 

/^-rix'^+2izx'^^f 

2i dz 


1 


^-rix''^+2izx' 


2ri 


= -- 2zF{z/^) + , (B.9) 


= — IZ, -e 


where we used dF/dz = —2zF -|- 1. 
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B.l Harmonic oscillator 


In the case of the harmonic oscillator on the half-line, the Feynman kernel is given by 

'imio{{x‘^ + x'‘^) cosut — 2xx'}' 


K{x, t; x', 0) = 


mu 


exp 


27rihsmut { \ 2h sin ut 

imu{{x‘^ + x'‘^) cos ut -|- 2xx'} 


=F exp 


2hsin ut 


(B.IO) 


where the minus (plus) sign corresponds to the Dirichlet (Neumann) boundary condition. 
From our initial profile with the Dirichlet boundary condition; 




the time-dependent solution is 




mu 


2r] V 27rihsmut 
1 


exp 


Xn imux^ COS ut 


4(7^ 


+ 


2hsinut 


X 




where complex functions t], and ( are dehned by 

1 imu cos ut 


7]{t) 


4(7^ 


,N ^0 -Po 
C{x,t) ■ = —+t — 


2h sin ut ’ 
imux 


2a2 


h h sin ut ’ 


^0 , -Po , imux 
C{x,t) :=— + t— + 


2a2 


h h sin ut 


K 




(B.ll) 


(B.12) 


(B.13) 


On the other hand, from our initial prohle with the Neumann boundary condition; 

'' {x-xof .poX 


$(a;,0) ={x - {^ + 1^ 


exp 


4(7^ 


+ F- 


h 


the time-dependent solution is 

T-—lexp 


X 


27iih sin ut 
I (1 pK 


Xn imux"^ cos ut 


+ 


[2p\2y p 

Xq , .Po 


4(7^ ' 2hsinut 

+ Ce^"/^A+tCF(- 


V 2Vh 


+iCF 


K 




+2 


1 pK 


2\1 p 


eCV4r) ^ gCV4^ ] + 




Vh V Vp \ ‘^Vp 


(B.14) 


iC 


(B.15) 
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B.2 Free particle 


The expressions for a free particle {y{x) = 0) are obtained in the limit a; —)■ 0 from the 
ones in the previons snbsection. The Feynman kernel is given by 


^ I ^im{x—x')^I2ht ^ im{x+x')^I2ht 

“T ^ 


^-ihtk?l2m, ^ ^ik{x-x') ^ ^ik{x+x') | 


1 


2tt 


(B.16) 


The time-dependent solntion from onr initial prohle with the Dirichlet bonndary condition 
is given as 


<F(x, t) = 


m 


^—aiQ / 4(T^ / 2ht 


2r] V 2Tiiht 
1 


X 






2^ 


(B.17) 


where complex fnnctions f], (, and ( are now 

... 1 im xo .po imx - Xq .Po , imx 

while the time-dependent solntion with the Nenmann bonndary condition is 


<h(x, t) = 


^ ^—x^jAcP'^ivnx^lXht 


X 


2niht 

2p \ 2 Y V 

Xo , .Po 


(£!«"/-” + Ce 


C^/4r? 






2x/h 


h V J Vv \ 2^y ^ 


V 2^ 


+- 


+2 

F 




2^/h 


(B.19) 
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